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The recently BaBar results raise the possibility that X(3872) has negative parity. This makes people recon- 
sider assigning X(3S12) to the ID2 (cc) state. In this paper we give a general form of the wave function of 
2 h mesons. By solving the instantaneous Bethe-Salpeter equation, we get the mass spectrum and correspond- 
ing wave functions. We calculate electromagnetic decay widths of the first 2 h state which we assume to be 
the X(3872) particle. The results are r(2~+(3872) -> J /\\iy) = 1.59^2 keV, r(2~+(3872) -)• \|f(2S)y) = 

2.87+J;^ eV and r(2~+(3872) -> \|/(3770)y) = OA35tom7 keV - The mtio of branch fractions of the second 
and first channel is about 0.002, which is inconsistent with the experimental value 3.4 ± 1.4. So X(3872) is un- 
likely to be a 2 h charmonium state. In addition, we obtain a relatively large decay width for 2 (3872) — > h c y 
channel which is k e V. 

PACS numbers: 13.25.Ft, 13.40.Hq, 12.39.Ki, 12.39.Pn 



I. INTRODUCTION 

In recent ten years B-factories have found many new heavy resonances which cannot be temporarily assigned to any char- 
monium or bottomonium states predicted by quark potential models. These particles are named as X, Y, Z particles. Among 
them X(3872) is the most famous one and the one which has been studied most carefully by experiment and theory. It is also 
the only one that has been detected through many decay channels. Since X(3872) was found by the Belle Collaboration jjj], 
other groups like the CDF 12!], D0 [3] and BaBar [4D all confirmed its existence. The following experiments put some effort into 
studying its quantum numbers. By detecting the EM decay channel X — > J/yj, the Belle group ]5[] fixed the charge parity of this 
particle to be positive. But the parity is disputable. The Belle's results J611 favor positive parity, while the CDF H concludes that 
both positive and negative parity are possible. Recently, by analyzing the channel X — > J/\\i(£> — > J /\\fK + K~n° with the entire 
BaBar data sample collected at the T(45') resonance, the BaBar Collaboration [8] concludes the negative parity for X meson is 
preferred. But just as Ref. [9] pointed out, positive parity assignment cannot be ruled out completely by the analysis. 

For the mass of X(3872) happens to lie around the DD* threshold, many authors believe the molecule state is found. Most 
of the work about this particle is built on the molecule assumption or the extension of this model to mix some cc or other 



19J, virtual state IBOl l2lh . etc 



components fl 101 - 41 €41 . Other fashionable models include hybrid state 11711 . tetraquark state 1118 , 
Although it's not as exciting as the former models, the traditional charmonium interpretation still needs careful investigation. 
This model faces two main problems. First, quark potential models have not predicted any state which has mass near 3872 
MeV. Second, the detected strong decay channels show there is large isospin violation. But these cannot exclude this assignment 



completely. As pointed out by Ref. 12211 . large isospin violation can be explained by the kinematic suppression effects, and the 
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mass of a charmonium may be affected if one considers coupled channel effects. 

Before the BaBar Collaboration gave their results H, most charmonium models assigned %ci (2P) to X(3872). Now negative 
parity is favored by the latest experimental results, so Dj state has to be reconsidered 112 311 . This state is the only ground state 
of spin-singlet D wave charmonia which have not yet been found. So finding this particle and determining its properties will be 



helpful to get information about interactions of quarks inside mesons. In Ref. |24j] , the authors have studied the 2 + production 



in semileptonic B decays. Because EM decay channels are clean, they are good channels to find new particles. It is important 



to calculate EM decays of 2~ + in order to fix down whether or not X(3872) is this state. Ref. 112511 used light-front model and 
Ref. 1 26] used pNRQCD to do the radiative transition calculations. Both of them got a negative result. Ref. 1270 studied the 
radiative transitions and the K°D°D decay mode. Ref. 12811 calculated the 2~ + (cc) production cross section at the CDF using 



fragmentation functions. All of them got results contradict the experimental value of X (3872) production. 



Beside the charmonium interpretation, Ref. fl29ll considered the 2 tetraquark model. Their conclusion is X(3872) cannot 



have a 2 tetraquark structure. Using heavy hadron chiral perturbation theory, Ref. 13011 discussed the molecule interpretation 



of X (3872) both in the 1 ++ and 2 + cases. But just as [22] pointed out that even in the 1 ++ case this one pion exchange bound 



state is dubious, so the molecule state with larger angular momentum will be strongly disfavored. 

In our previous paper l 3lL we have investigated the case which X(3872) is assumed to be the % c \ (2P) state. The EM decay 
ratio we get consists with the BaBar value 13211 . There we also get Br(% c i (IP) — >J/\^y) = 35.6% (this value changes to 42.4% by 
setting the parameters equal to new values in this paper) which is close to the experimental value (34.4 ± 1 .5)%. We also studied 
the J/y— > TjcY process to make sure the method is reasonable. The branch ratio is 2.2%, while the PDG value is (1.7 ±0.4)%. In 



this paper we will use the instantaneous Bethe-Salpeter (BS) (33 



3511 method to study the 2 + scenario for X (3872). Although 



other people's work disfavors this assignment, we think it's still necessary to give it a careful study with a different method. First, 
for the radiative channels different models got very different results. On one hand the discrepancy comes from the differences 
of these models; on the other hand people usually concentrate on determining the quantum number of this particle, thus rough 
estimations tend to be used. But in the radiative transition processes, the gauge invariance which is important may be violated 
to a large extent by using some approximations. Therefore, a careful study is needed. Second, even X(3872) is not l l D 2 (cc), 
our study will provide some useful information for this undiscovered charmonium. This will be even helpful if the 2 state 
has similar mass with that of 2 h which will bring more challenge to the experimental detection. Quark potential models have 
predicted that the mass range of 1 '/^(cc) state is 3760 ~ 3840 MeV [26] which is below the DD* threshold. This will lead the 
particle to have a narrow decay width. Here we will not consider the threshold effects but adjust the parameters to fix the mass 
of ground state of 2 h around 3872 MeV. 

This paper is organized as follows. In the second part we give the general form of the wave function of 2~ + states. Then we 
present the instantaneous BS equation which satisfied by the 2 h wave function. By solving the coupled equations we get the 
eigenvalues and corresponding wave functions. In the third section with Mandelstam formalism we calculate electromagnetic 
decay widths of this state by assuming it has the same mass with X(3872). In the fourth section we give our discussions and 
conclusions. Appendix shows the positive wave functions and details of form factors. 
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II. INSTANTANEOUS BS EQUATION OF THE WAVE FUNCTION OF 2 + STATE 



The wave function of 2 + states with mass M, momentum P and polarization tensor e^v has the general form 



M 



M ' 



(1) 



which satisfies constraint conditions of the Salpeter equation 113411 . and then we obtain the following relations 



Mq±) = 



Mq±) 



fi(q±)M(mi(02-m2(Qi) 
-f2(q±)M((Oi+(0 2 ) 



(2) 



(miC02 + C0l«l2) 

where q is the relative momentum between constituent quark and antiquark which have masses m\ and m 2 , respectively. q± is 



defined as q — jpP and ft),- has the form y mf — q 2 ± . f\ ~ f\ are functions of \q\. One can see when m\ = m 2 , the particle has 
definite C-parity. The term with fj which has negative C-parity disappears with the equal mass condition. So the wave function 
just has two independent quantities, f\ and f 2 . 



With the same method used in Ref. 13611 we can get the coupled instantaneous BS equations for 2 + state which have the 
following form 

CO] 



(M-2& l )(f l (q) + —f 2 (q)) 
mi 

M 3 - 1 



tf.kf - -q 2 P] [ mi (V s - V v )( mi f 2 (k) + COj/i (k)) - (V s + V v )k ■ qf 2 (k)} , 
lq*m\(S)\ i 

ft>i 



(3) 



{M + 2m)(fl(q)-—f2(q)) 
mi 



= - / d'k 



^^—[{q-kf - \f&\ [mi (V s - V v )(mif 2 (k) - mfi (%)) - (V s + V v )k ■ qf 2 (k)} . 

To solve above equations, we have used the Cornell potential ( This phenomenological potential already reflects the main 
feature of the interaction between quark and anti-quark. One can modify this potential by introducing additional high order 
terms to get better spectra, but it will be hard to solve the BS equation. Because the wave function constructed is the most 
general one, even by using this simple potential we can get a reasonable result ) which in momentum space can be written as 

V(q) =V s (q)+y ®fV v (q), 



V v (q) = - 



2 a s (q) 



(4) 



3n 2 q 2 + a 2 '' 
12k 1 



27 



ln(a + /— ) 



QCD ' 



Here we adopt the following values for the parameters |l37|, a = e = 2.7183, a = 0.06 GeV. X = 0.21 GeV 2 , m c = 1.62 GeV, 

flDB, 



Aqcd = 0.27 GeV. One notices that we used different values of these parameters as that in 



for these new values can 



adopt to more mesons with different quantum numbers 113 711 . even though the former can lead to better spectra. m c is the mass 



for the constituent charm quark, which is a little larger than half of the mass of r| c or caused by the small binding energy. 



Because we want to study electromagnetic decay of 2 h (l£>) as X(3872), we adjust Vo to make the first state with this quantum 
number to have mass equal to 3872 MeV (Since we do not consider higher order interactions, to get the mass splitting between 
mesons with different quantum number we introduce Vq to fit data. Changing other parameters will affect the spectrum a lot, 
while changing Vo just cause a translation of the spectrum). By dong so, we get Vo = -0.044 GeV. 
The normalization condition for the BS wave function is 

f _4_Tr[9 ++ ^( P ++ --r--'p--]=2M, (5) 
J (2%) 3 M M M M 

where cp ++ and <p are positive and negative energy part of the wave function, respectively. Their expressions can be found in 
the Appendix. Putting wave functions of 2 h into above equation, we get its normalization condition: 

f d 3 q 8 COi M 

To solve Eq. (fjjl, we first discretize the relative momentum q and k. For the wave function will approach when \q\ becomes 
large, we cut off \q\ at 7.71 GeV (where the wave function is small enough). By solving the eigenvalue equation we get the mass 
spectrum and corresponding wave functions (in our results there is no degeneracy). Masses of the leading four states have been 
listed in Table|T] 

III. EM DECAY OF X(3872) AS THE l 1 D 2 {cc) STATE 



The wave function of the 1 state has the form 113 811 



<Pi- (9x) - (q± -e)[/i (q±) + jjMlx) + + -jgrMlx)] +Mt]f 5 (qx) + ^/e(?±)] 



(7) 



where £ is the polarization vector of the meson, and /,s are scalar functions of q. When m\=mj_, f% and fi equal to 0, which 
makes the wave function have negative C-parity. 

The electromagnetic transition (see Fig. 1) amplitude is 

T = (P f e 2 ,ke\S\Pei) = ^ d\P f +k-P)e^, (8) 

where e q = I is the charge of the charm quark in units of e; e, £i and £2 are the polarization vectors (tensor) of the photon, the 
initial meson and the final meson, respectively; 94^ is the hadronic transition matrix element. With the method which has been 



proved to be gauge invariant in Ref. II39I1 . at leading order can be written as (Here we only keep the terms contain positive 



energy wave functions. Other terms only contribute less than 1%.) 

M% = 1 7^ Tr [^' ++ (^ + a 2^±)^9 ++ (?±)-9' ++ (?±-ai^)^ ++ (^)Y 4 ], (9) 

where Pfj_ = Pf— -jgfP> 9 ++ an d <p /++ are the positive energy wave functions of initial and final particle, respectively. (f> ++ is 
defined as 7°(cp ++ )^7°. In the charmonium case, OCi and 0C2 equal to i. 



After doing the trace and integrating out q, we get the following form of the amplitude 

9vfi = e a p£ M P a P f5 (e& a5 P?pfti + e^ oSa P^f 2 + £^M 2 g a ^t 3 ). (10) 

One can see that there are three form factors t\ ~ t 3 which are integrals of q. Their explicit expression can be found in the 
Appendix. We give the results of these form factors for different states in Table [TT] Eq. ( TTOb has a different form with that in 



Ref. 12511 . but one can check they are actually equivalent. 



Another possible EM decay channel is 2 + — > 1 + y. Similar to the above calculation, we first present the wave function of 



1 + " |40[] 



q>i+- (ax) = {qx-m («x) + jjfi(q±) + ^ftiix) + -^tMix)]y 5 . (U) 

One notices that it has the same Lorentz structure inside square brackets with the wave function of 2 '". From Eq. <j9j we get the 
expression for the amplitude of this transition 

= e^M 2 (P^P}s l + P^ 2 ) + e^M 4 s 3 (12) 

This amplitude which contains seven form factors is a little bit complex compared with Eq. (TTOb . The procedure to calculate 
these form factors is similar to that of above, we will not give their explicit expressions but just list their values in Table Hill 

TABLE I: Masses (in unit of GeV) of l D 2 {cc) mesons with Vo = -0.044 GeV. The uncertainties are gotten by varying all parameters in Cornell 
potential by 5%. 



n l D 2 (cc) 


ID 


2D 


3D 


4D 


Mass 


o 079+O.179 
J - o/z -0.180 


4 , 98 +0.189 


4 450 +0.196 


4 655+ 201 



TABLE II: Form factors (in unit of GeV 2 ) of the decay process 2 + (3872) — ¥ \\t(nS)y. The uncertainties are gotten by varying all parameters 
in Cornell potential by 5%. 

form factor t\ t 2 ?3 

2-+ (3872) ->//Vy 
2 _+ (3872)-H|f(2S)y 
2-+ (3872) ->\|/(3770)y 



8 36 +0J0 447+ - 047 -4 37+0-58 , Q -2 

Oi:,D -0.62 Ui ^'-0.040 ^""-0.68 x lu 

23.5+^f, - 4 - 47 ^L4o x 10-2 3 - 34 -a4i x 10-2 

7.3+1.97 ]4 9+1.7 _0 9 98+ 084 

'■■"_0.64 ^- z -1.2 u - w6 -0.096 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

By solving corresponding BS equations we also get the wave functions of 1 states and l + ~ states. We use the same 
parameters as that used when we solve the equation of 2 h except Vo. By fitting the mass spectra, we find the best-fit values of 
V : -0.144 GeV for 1+" and -0.1756 GeV for 1~ . 
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TABLE III: Form factors (GeV 2 ) of the decay process 2 +(3872) — > h c (lP)y. The uncertainties are gotten by varying all parameters in 
Cornell potential 5%. 



form factor 




S2 


S3 


i 4 


■55 


i'6 


*7 


2-+ (3872) ->h c (lP)y 


-3 96+ 50 


J - z -0.71 


-0 742+°- 117 
u - '^ z -0.145 


2 44+ 038 
z -^ -0.31 


20.81^ 


28.6^1-9 


4 4Q+0.73 




m 1 



P> 



Pi 



FIG. 1 : Feynman diagrams of the EM decay of 2 + meson. 




Using Eq. ( fTOb we can calculate decay widths of 1 7>2 — > \y(nS)y, which have been listed in Table JV] We also calculate 
the uncertainties caused by varying the central values of the parameters in the Cornell potential by 5%. The form factors are 
presented in Table [TT] 

Considering the upper limit of the decay width, r X n872) < 2.3 MeV, we can get the theoretical predictions for the lower limit 
of branch fractions of these channels: Bi(X ->■ J/yy) > 6.87 x lfr 4 , Bi(X -> w(2S)y) > 1 .25 x 1(T 9 , Br(X -> v|/(3770)y) > 



5.87 x 10" 5 . There are experimental lower limits on the first two decay channels |41j, which is much larger than our predictions. 
One can see that the branch fraction of X to J/\\t is one order of magnitude larger than that of X to \|/(3770), while it is much 



(13) 



larger than that of X to \|/(2S). This prediction is strongly in contradiction to the experimental value 113211 

Br(X(3872)^ W (2S)) 
Br(X(3872)^y//» 

To understand our results, in Fig. 2, we plot the wave functions f\ of 2~ + state, of \|/(3770), /5s of //ty and \|/(25) (For 2 
state, /1 almost equals to fa. For J/\\i and \|/(25), /s and fa which have almost equal absolute values give the main contribution, 
while for V|l(3770) fa and fa are the large part of the wave function.). To make the wave functions dimensionless we multiply 
them by different quantities. One can see that the wave functions of 2 ^ reach their maximum values when |^| takes about 0.75 
GeV. Also, the wave function of \|/(3770) gets its maximum value at the same \q\. This state which is the third eigenstate of our 
1 instantaneous BS equation includes 5-wave and D-wave mixing. For this state, one can see q 2 fa has the shape of D-wave. 

Actually our results above can be easily understood from the node structure of wave functions of initial and final particles. 
One can see there is a node in the wave function of \|/(25), which locates almost at the same position of \q\ where the wave 
functions of 2~ + gets their maxima. This causes strongly decreasing of the overlap integral. While for J/\\t and V|l(3770), there 
is no node, so we get a much larger width. We have used the same arguments in our previous work to explain the results of 1 ++ 
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-20 




q *f3of\|/(3770) 
M 2 *f5of\]/(2S) 
M 2 *f5ofJAt/ 
M*q 2 *fl of 2" + (3872) 



0.5 1 1.5 2 2.5 3 

q (GeV) 

FIG. 2: Wave functions of 2~+(3872) and of J /\f, \|/(25), \|/(3770). 



3.5 



decay [3JL] . Because the wave function of 2P state has a node, the cancellation happens for 15, not for 25, which leads to a result 



closer to experimental value. 



One can see the decay width of the J/\\tJ channel is close to that with other methods, while the \|/(25)y channel has a quite 



different width. Ref. 125 



26 



14311 consider S-D mixing by introducing a mixing angle, say 8 = 12° in Ref H25LI26H . Even with 
the same mixing angle there is still a large discrepancy which manifests the sensitivity of the results to the node structure. As for 
our method, all the adjustable parameters are contained in the potential and the wave function are constructed in a general way 
which makes it contain mixture automatically. 

If we trust the results of BaBar H, there will be no candidate for this particle in the spectrum predicted by quark potential 
models. The inconsistent of theoretical predictions and experimental data will force us to abandon the 1 '/^(cc) assignment to 



X(3872). But we think it's still too early to abandon the charmonium assignment. Just as Ref. 1421 pointed out that experimental 
data interpretation is a subtle issue. We have to make sure what exactly we have measured, e.g. is there any possible that two 
particles with masses close to each other have been detected, such as 2~ + and 2 ? 

If l l D2 (cc) proved not to be a candidate of X(3872), we would like to give predictions of this missing state. By setting Vo 
= -0.113 GeV, we resolve the coupled equations satisfied by the 2~ + state. Then we get the ground state with M= 3820 MeV 



and corresponding wave functions (Quark potential models give the mass range 3760^3840 MeV Ref. 12611 . We choose 3820 
MeV as an example). The decay widths of the particle with this mass are listed in Table II VI One notices that T(1 1 D2 — ^ J/WY) 
changes only a little, while for the other two channels the widths change 4^10 times. This comes from the sensitivity to the 
node structure and phase space changing for the last two channels. 

For the decay channel 2~ + (3872) — > h c (lP)y, we also use two groups of mass values to calculate the width. For M=3872 
MeV, the decay width is 391 keV. Within the uncertainty range, our result is consistent with those from other approaches (except 
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Ref. 12611 where the relativistic corrections are neglected). ForM=3820MeV, the result is 395 keV which just has a little change. 



This is different to that in Ref. 14511 . which changes more than 100 MeV when the mass of the particle changes. We can see 
this channel is more than two orders of magnitude wider than 7/\|/y channel, which means this is the main EM decay channel of 
l l D z . 

In summary, we have gotten the mass spectrum and wave functions of 1 Z?2(cc) states by solving the corresponding in- 
stantaneous BS equations. Then within Mandelstam formalism we calculated decay widths of different EM transition pro- 
cesses: r(2-+(3872) -> 7/\|/y)=1.59 keV, r(2~+(3872) — » \|/(2%) = 2.87 eV, r(2-+(3872) -> \|/(3770)y) = 0.135 keV and 
r(2~ + (3872) — > h c y) = 392 keV. The ratio of branching fractions of the first two channels shows that X(3872) is not likely to 
be a 2~+ state. More precise measurements and analyses are definitely needed for the decay widths to clarify the nature of this 
mysterious particle. 



TABLE IV: Electromagnetic decay widths of l l D2{cc). For our work, the results out (in) the parentheses are gotten by setting M x = 3872 
MeV (3820 MeV). Ref. H] used M x = 3872 MeV (out parentheses) and M x = 3837 MeV (in parentheses). Ref. j^l and Ref. |H] did 
the calculation by setting M z =3872 MeV, while Ref. (H, Ref. (H and Ref. 03] did the calculation with M x = 3820, 3825 and 3796 MeV, 
respectively. The uncertainties are gotten by varying all parameters in Cornell potential 5%. 



Ref. 



r( 1 Z) 2 ^/vY)(keV) T( 1 D 2 ^J/WY) (keV) 



r( 1 £» 2 ^ V (25)Y)(eV) 



r(»» 2 -M|/(3770)y) (keV) 



This work 
Ke&Li [25] 
JiaetaJ [26] 
S&Z[43J 
Eichten et al [44J 
B&G[45J 
Chao et al [46J 



392 



+62 
411 



(395) 



587 ~ 786 

288 

303 
464(344) 

375 



1.59^3(1.08) 
3.54 
3. 1 1~ 4.78 
0.699 



2.87+^(0.682) 
0.60 
17—29 
1 



0.135+^(0.0128) 
0.356 
0.49-0.56 

0.34 
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APPENDIX 



The definition of (p++ can be found in Ref. 13611 . From the general form of wave functions we get the positive energy part of 
wave functions of 2 h , 1 and 1^ : 



cp++ + = e^q v ± (Aij 5 +AoJj 5 +A 3 fy y 5 ), 



(14) 



q>+t =Btt+B 2 tf + B 3 ?fy[ x + B 4 q± ■ e + B 5 q ± ■ zf + B 6q± _ ■ e^ ± + B n q L ■ e} 



(p+t = q± -e(CiY5 + C 2 H +C 3 ^ ± Y 5 ), 



(15) 
(16) 



where A,s, B,s and C;S are scalar functions of q which have the following forms 

Ai = -(/i + — / 2 ), A 2 = — — / + — f 2 , A 3 = -— — + — / 2 ), (17) 
2 mi 2Mcoi mi 2Mcc>i mi 

d M f coi mi CO! 1 coi 

0\=—(j5 /6j, B 2 = - — {f5 J6), a 3 = —[f 5 J6), 

2 mi 2C0i mi 2C0i mi 

Mf mi if 2 /ni . 1 (Oi 

^4 = ^/5-7-/6-^ (/3 + — A), fi 5 = --— /5 /6 , (18) 

2mi C0i lMfm\ C0i 2(0i mi v 7 

D M / , 1 it , "Im D 1 , , m i /, , m i r \ 

B6 = " 2mT^ /6 + 2M^ (/3 + W /4) ' B? = " ^W /5 + 2l^m7 (/3 + ^ 

d = i(/ 1 + % 2 ), C 2 = -^-(/ 1 + ^/ 2)l C 3 = --!-(/! + ^/ 2 ). (19) 
2 mi 2Mcoi wzj 2Mcoi mi 

To get the expressions for <p , one just need to change C0i to — C0i in cp ++ . 
The form factors for the 2 h — > 1 y process are expressed as 

h = (ai+b 1 )-(a 3 +b 3 )- a 2 P ■ P f (a s + b % ) + a 2 P- P f (a 12 + b n )-2P- P f (a 15 - a i6 ), 

P ■ Pf P Pt 

t 2 = -2(a 2 + b 2 ) + 2a 2 — ^-(a 5 -b s ) + 2a 2 P-P f (a n +b n )- 2{a 2 —^-) 2 {a u -b u )i (20) 

2 , , P Pf , 

h = -j^{ai4 + bu)+2a 2 — r (an+bn), 



where a,s and bjS are the integrals of q: 



d 4 A T3 i A D ^ l^l 2 

f 

4 



/ (2^3 AljB2(A2jB3) 2^ (3C ° s2e_1) ' 



^ (2% AijB7(A3jB7) 8S|^P (5 



fl3(«i2)= / 7 ^AiB 7 (A 3 B 7 )--^ T (5cos 4 e-6cos 2 e+l), 
J \2%y 8|P/| 2 



«4(ai3)= / -^ T AiB 7 (A 3 B 7 )^(cos 4 e-2cos 2 e+l), (21) 
J \2Tiy o 

, ^ f dq , . lol 3 , . 

05(07,010,014) = / 7r-T T AiB 7 (A 2 B3,A3B 2 ,A 3 B 7 )^^-(cos^ 9-cose), 
J (2ny 2\Pf\ 

a 6 (a 9 ) = -J ^L A2 B 3 (A 3 B 2 ) ^^j(5 cos 3 9 - 3cos9), 

«i5(ai6)= / 7-^-3A 2 B3(A 3 B 2 )-^-(5cos 3 e-3cose), 
J (2ny 2|P/| 3 

where is the angle between Pf and q; the arguments of A,- and B, are g and q + a 2 Pf, respectively. To calculate b u one just 
need to change the argument of B, from q + a 2 Pf to q — OL\Pf. 
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